Viscoelastic Behavior of Solid 4 He 
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Over the last five years several experimental groups have reported anomalies in the temperature 
dependence of the period and amplitude of a torsional oscillator containing solid 4 He. We model 
these experiments by assuming that 4 He is a viscoelastic solid — a solid with frequency dependent 
internal friction. We find that while our model can provide a quantitative account of the dissipation 
observed in the torsional oscillator experiments, it only accounts for about 10% of the observed 
period shift, leaving open the possibility that the remaining period shift is due to the onset of 
superfluidity in the sample. 
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In 2004 Kim and Chan pi Q reported anomalies in the 
resonant period of a torsional oscillator (TO) containing 
solid 4 He. With exquisite sensitivity they detected a re- 
producible decrease in the oscillator period upon lowering 
the temperature below 200 mK. Subsequent experiments 
in several laboratories [1 S, I, 0, S S, S EH EH have now 



reproduced these results. The size of the effect depends 
on the sample preparation, pressure, and concentration 
of 3 He impurities. References pj| and p]| review the 
recent theoretical and experimental work. 

A natural interpretation of the TO anomalies is the on- 
set of the elusive and long-anticipated supersolid phase 
of matter [3, In a supersolid, superfluidity co- 

exists with the crystalline order of a solid; one expects 
a supersolid to exhibit superfiow, much like a supcrfluid. 
Leggett pj| proposed that the superfiow is best detected 
by searching for "nonclassical rotational inertia" ; a su- 
perfluid condensate would remain at rest and not par- 
ticipate in rotation, and the resulting mass decoupling 
would reduce the rotational inertia and decrease the res- 
onant period of oscillation. While compelling, the su- 
persolid interpretation of these experiments has yet to 
be corroborated by other measurements, such as the re- 
sponse to pressure differences Ijj. Moreover, Day and 
Beamish |18| reported a pronounced increase in the shear 
modulus of 4 He at temperatures below 200 mK, with a 
dependence on measurement amplitude and 3 He impu- 
rity concentration similar to the TO anomalies. Their 
results suggest that changes in the shear modulus might 
be intimately related to the TO anomalies. 

This Letter presents a detailed discussion of the me- 
chanical response of a TO containing a viscoelastic solid. 
We build on earlier work by Nussinov et al. [19( , who cor- 
rectly identified a "back action" term in the TO equation 
of motion that represents the dynamical effect of the solid 
helium on the torsion cell. However, in contrast to Nussi- 
nov et al. we find no need to assume that the solid helium 
behaves as a glass. Instead, with a few carefully stated 
assumptions, we find that we can model the solid helium 
as a classical viscoelastic solid — i.e., a solid with internal 
friction. The TO period shift and dissipation peak are 



naturally related to the real and imaginary parts of the 
frequency-dependent shear modulus of the solid. We use 
our results to fit the dissipation peak reported in Clark et 
al. [Bj], and extract a temperature dependent time scale 
t(T) from the data. With all of the phenomenological 
parameters determined, we find that we are only able to 
account for about 10% of the period shift reported in 
Ref. [H, leaving open the possibility that the remaining 
shift is due to the onset of superfluidity (or supersolidity) 
in these samples. 

Following Nussinov et al. [13], we assume that the 
empty torsion cell is perfectly rigid, with a moment of 
inertia J sc about its rotation axis. For small angular dis- 
placements 9 the torsion rod provides a restoring torque 
— a9, with a the torsional spring constant. There is also 
a damping torque — r y osc d9/dt, with 7 osc a temperature- 
dependent dissipation coefficient. The cell is driven by 
an external driving torque of T ex t(f)- Finally, the solid 
helium inside the torsion cell exerts a moment M(t) on 
the cell. The equation of motion for the cell is then 
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t (f) + M(i), (1) 



with M(t ) — J dt'g(t — t')9(t') for a linear system that is 
invariant under time translations. We can Fourier trans- 
form the equation of motion to find the response function 
x(Lu)=9(uj)/r cxt (Lu) of the TO: 
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(2) 



The response function is of fundamental importance in 
interpreting the TO experiments, as its poles give the 
resonant frequencies and their quality factors. 

The complex back action term g(ui) contains all of the 
information about the dynamical response of the solid 
helium, and modeling this quantity is the focus of the 
remainder of the paper. Before delving into calculational 
details we can make a few general statements about g(to). 
First, causality requires g{z) to be analytic in the upper 
half of the complex z-plane, so the real and imaginary 
parts of g(z) satisfy Kramers-Kronig relations for z — lu. 
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Second, if the helium behaves as a perfectly rigid, nor- 
mal solid, g(uj) = In e cu 2 , with Jhc the rigid body mo- 
ment of inertia of the helium. As we will show below, a 
finite shear modulus produces corrections to this result 
which vanish with a higher power of u> at low frequencies. 
Therefore g{u> = 0) = 0; the back action only modifies 
dynamical quantities. At low frequencies the "glassy" 
model for g(u>) proposed by Nussinov et al. [see their 
Eq. (20)] goes to a constant and induces an unphysical 
shift in the static spring constant a. 

To calculate the back action term we first need to 
model the properties of the helium in the torsion cell. 
To keep the description as general (and simple) as pos- 
sible we model the helium as an elastic continuum, with 
an equation of motion 



pc%Ui = djaij, (3) 

where p is the density, itj is the i th component of the 
displacement field and is the stress tensor. The stress 
tensor has both a reversible piece and a dissipative 



piece afj, with ct„ _ „ tJ , „ y 



+ a ij |20j ■ For a linear medium 
the reversible stress is linearly related to the strain tensor 
"a- = {d k Ui + dtUk) /2 by = \ i3 kiu k i, where X i3k i is a 
fourth-rank tensor of elastic coefficients. For a uniaxial 
crystal such as the hep phase of 4 He there are 5 inde- 
pendent components of Xijkf, however, to simplify the 
analysis we assume that the helium can be modeled as 
an isotropic elastic solid, with cr^ = 2/zou.y + XoSijUkk, 
where Ao and po are the Lame coefficients (with po the 
shear modulus). The dissipative part of the stress ten- 
sor, which describes the internal friction of the solid, 
must be odd under time reversal, and can only depend 
on gradients of the velocity Vi — d t Ui [13]. For a linear 
medium af 3 = rjijkiVki, where rjijM is the viscosity tensor 
of the solid and Vik = (diVk + dkVi)/2. Again assuming 
anisotropic medium, afj — 2rjVij + (( — 2ri/3)5ijVkk, with 
r) and £ the shear and bulk viscosities of the solid. With 
these simplifications the Fourier transformed equation of 
motion becomes 



- puj 2 u = B(uj)V(V ■ u) - p(uj)V x V x u, 



(4) 



with B(lo) — Xq + 2po — iuj{C, + 4/7/3) and p{uj) — 
po — iurj = po(l — iTU)), with timescale r = i]/pq. In pass- 
ing we note that this model, known as the Kelvin- Voigt 
model, is among the simplest of viscoelastic models — a 
single "spring" (the elasticity) is in parallel with a sin- 
gle "dashpot" (the viscosity). More elaborate models, 
involving series and parallel combinations of springs and 
dashpots, can produce a shear modulus p{u) with a more 
complicated frequency dependence. For an example of a 
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similar analysis for colloidal crystals, see Rcf. 

The next step is to determine the response of the he- 
lium inside the torsion cell to the rotation of the cell. 
For simplicity we will present results for a torsion cell 
that is an infinitely long cylinder of radius i?; the gener- 
alizations to an annular geometry or a cylinder of height 



h are straightforward 22J. If we assume the helium is 
in perfect contact with the walls of the torsion cell (no- 
slip boundary conditions), and the torsion cell oscillates 
about its azimuthal axis with a frequency to and ampli- 
tude 9$, then the induced displacements in the helium 
have the form u = ug(r)e~ luJt 9. Substituting this into 
Eq. (|4]) and solving the differential equation with the no- 
slip boundary condition ug(r — R) = R9o, the solution 
that is finite at r = is 



u e {r) = R9 Q 



Ji(kr) 
Ji(fci?)' 



(5) 



where k 2 — uj 2 pj p(uS) and J\{z) is the Bessel function of 
order 1. In this geometry the torsion cell only induces 
shearing displacements in the helium. 

The final step of the calculation is to determine the 
moment that the oscillating helium exerts back on the 
torsion cell. The only non-vanishing component of the 
stress tensor is oe r — p{u){d r — l/r)ue(r); evaluating this 
on the surface of the cylinder, integrating over the area 
of the cylinder to obtain a force, and then multiplying by 
the radius to obtain a torque, we find the moment 



M(t) 



a U 2 (kR) 
-t) uj -f H c e 



kRJx{kR) 



(6) 



where Thc = irphR i /2 is the rigid body moment of inertia 
for the helium inside the torsion cell. In terms of the 
back-action term g(ui) defined in Eq. ([2]), 



g(oj) = Ih c u> 2 + In c uJ 2 g(kR), 



where the function 



4J 2 (z) 
xJ\{x) 



1. 



(7) 



(8) 



is the correction to the rigid body result due to the finite 
shear modulus of the helium. 

To simplify our result, we note that for a typical TO 
the speed of transverse sound ct = ypo/p ~ 270 m/s, 
the frequency / ~ 10 3 s -1 , and the radius R ~ 0.5 cm, 
so that \k\R = 2irfR/cT ~ 0.1. Therefore we can safely 
expand Eq. JHJ for small x, with the result g(x) ~ x 2 /24; 
then Eq. ([2|) becomes 



-it, 



t U! -i^f osc uj + a- 



pR 2 Lu i I Hc F(R/h) 
2Ap(oj) 



. (9) 



where /tot = ^Ho + ^osc- We see that the last term in 
Eq. ([9]) is the correction due to a finite shear modulus; 
for a perfectly rigid body, p — *■ oo and this term vanishes. 
Also, in this last term we have introduced a function F{x) 
to describe the effect of a finite cylinder height h [22| ; this 
function only depends on the aspect ratio x — R/h, with 
the explicit form 

192 °° 1 / \ 

^) = -^E(^i F < i ( 2 -- 1 H- (10) 
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For the infinite cylinder F(0) — 1, and more generally 
< F(x) < 1. For large x, F(x) ~ 2/x 2 - 744C(5)/7r 5 x 3 
with (,{n) being Riemann-Zeta function. In the particular 
case h = R, F(l) = 0.527. 

We now examine the effect of the viscoelasticity of 
the helium on the period and Q-factor of the oscillator 
by finding the poles of the response function, Eq. ([9]) 
(our analysis is similar to the procedure performed in 
Nussinov et al. [H|). Since = luqR/ct ~ 0.1 and 
Itie/hot ~ 10 -3 , we can treat the fourth term in Eq. (J9j) 
as perturbation about the rigid-body behavior, which 
has a resonant frequency loq — \J a/ /tot and dissipation 
Qq 1 = 7osc/\/^tot«- Expanding the poles about wo and 
Qq" 1 , and recalling that fi(uj) = /xq(1 — iro;), we obtain 
the fractional period shift AP/P = (P — P )/P and the 
shift in the dissipation AQ^ 1 = Q~ x — Qq 1 



AP 



A Re 



IM) 



AQ- 1 = 2 Aim 



Mo 
H(w ) 



A 



l + (ra;o) 2 ' 



= 2A 



1 + (™ ) 2 1 



(11) 



(12) 



where the dimensionless amplitude A is given by (recall 
ct = VW/ 9 ) 



.4 



48 /tot V ct 



(13) 



While v4 depends on the material parameters and the 
sample geometry, it is independent of the relaxation 
time t. For a typical TO the amplitude A is of order 
10~ 6 — 10~ 7 , so the resulting shifts are small, as as- 
sumed. As we will show below, amplitudes in this range 
can quantitatively fit the dissipation peak, but are a fac- 
tor of 10 too small to fit the period shift of the TO results. 

The simple Lorentzian form of these results suggests 
a strategy for fitting the TO experimental data. Notice 
that AQ -1 will have a peak as a function of tempera- 
ture T if the relaxation time r passes through the time 
scale loq 1 as the temperature is lowered; the peak occurs 
at T* such that loqt(T*) = 1, and at this temperature 
AQ _1 (T*) = A. Therefore, A can be directly determined 
from the peak value of the dissipation; we can then solve 
Eq. (fl"2"|) for luqt in terms of AQ~ 1 (T)/A as a function 
of temperature, allowing us to determine t(T). Having 
determined both A and t, we can then calculate the pe- 
riod shift AP/Pq as a function of temperature, with no 
additional fitting parameters. In passing we note that the 
temperature dependence of the shear modulus itself [HJ 
has a much smaller effect on the period shift and dissipa- 
tion than the temperature dependence of the relaxation 
time t. 

In Fig. [1] we show the relaxation time obtained from 
the measured Q-ia c to r of the BeCu TO used in Ref. [5| 
(to collapse all of the data onto a single curve we have 
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FIG. 1: (color online). Relaxation times derived from 
AQ _1 (T) of BeCu TO of Clark et al. @|. Green squares 
are for the blocked capillary (BC) sample, blue circles the 
annealed blocked capillary (ABC) sample, and pink triangles 
the constant temperature (CT) sample. The red line was 
found by fitting to all three samples and using t(T)luq = 
a exp(6T*/T)/[l + cexp(dT*/T)] with fitting parameters a = 
1.75 x 10" 3 , b = 7.55, c = 1.62 x 10" 3 and d = 7.01. 



scaled the temperature by the peak temperature T*). On 
this log plot we clearly see activated behavior both be- 
low and above T* , but with different activation energies. 
To account for this behavior, we will assume that the 
relaxation time r has the functional form 



r(T) 



r exp(E /T) 



(14) 



l + <5exp(£a/T) 

For BeCu TO data of Ref. @] we obtain t = 260 ns, 
6 = 1.62 x 10~ 3 , £ = 7.55T*, and E x = 7.01T*. At 
high temperatures, the activation energy for the blocked 
capillary sample of BeCu TO of Ref. [5j is found to be 
E = 260.4 mK (417.5 mK for the annealed blocked capil- 
lary sample and 341.4 mK for the constant temperature 
sample), and at low temperatures E = 18.6 mK (29.8 
mK and 24.4 mK, respectively). By using the derived 
relaxation time r the fits to the dissipation peak and pe- 
riod shift data from Ref. 0] are shown in Figs. [2] and [3] 
Having fit to the dissipation peak, we see that the same 
set of parameters accounts for only 10% to 20% of the 
period shift observed in these experiments. 

In summary, we investigated the viscoelastic behav- 
ior of solid 4 He at low temperatures. The response of 
a viscoelastic solid to an oscillatory shear stress is de- 
termined and used to study the anomalies in the reso- 
nant period and the dissipation of the TO experiments. 
Our approach is quite general, invoking the linear re- 
sponse of the helium together with the simplifying as- 
sumption of isotropy; our framework can be used to study 
other phcnomenological models for the mechanical be- 
havior of solid helium. For instance, we could also treat 
a model with a distribution of relaxation times, such as 
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FIG. 2: (color online). Peaks in the inverse of the quality 
factor of BeCu TO from Clark et al. 5]. Green squares 
are of the blocked capillary (BC) sample, blue circles the an- 
nealed blocked capillary (ABC) sample, and pink triangles 
the constant temperature (CT) sample. The red line is of the 
viscoelastic model with the change in Q-factor Eq. (|12[) and 
the derived relaxation time Eq. (|14|l . 
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FIG. 3: (color online). Resonant period shifts of BeCu TO 
from Clark et al. [|| . Green squares are of the blocked capil- 
lary (BC) sample, blue circles the annealed blocked capillary 
(ABC) sample, and pink triangles the constant temperature 
(CT) sample. The red line is of the viscoelastic model with 
the resonant period shift Eq. and the derived relaxation 
time Eq. (TJJ} . 



a "glass" model for the shear modulus [l9( of the form 
(i(u>) — /io(l — iruj) 13 . The simple single relaxation time 
Kelvin- Voigt model identifies a time scale associated with 
the viscosity of solid He; upon lowering the temperature, 
this relaxation time grows rapidly and eventually passes 
through ujq 1 , inducing changes in both the dissipation 



and the oscillator period. While the dissipation peak can 
be explained completely using the viscoelastic model, the 
model accounts for only a fraction of the period shift ob- 
served in Ref. [5} (although fits to some data [|| yield a 
smaller discrepancy between the model results and mea- 
sured period shifts). As originally suggested [l], 0], the 
remaining period shift may indeed be due to the onset of 
some type of superfluidity in the solid helium. 
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